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A 20 % substitution of Bi with La in the perovskite Bii-a^LajiFeo.sSco.sOs system obtained un¬ 
der high-pressure and high-temperature conditions has been found to induce an incommensurately 
modulated structural phase. The room temperature X-ray and neutron powder diffraction patterns 
of this phase were successfully refined using the 7mma(0,0, 7 )s 00 superspace group (7 = 0.534(3)) 
with the modulation applied to Bi/La- and oxygen displacements. The modulated structure is 
closely related to the prototype antiferroelectric structure of PbZrOs which can be considered as 
the lock-in variant of the latter with 7 = 0.5. Below 220 K, the neutron diffraction data 

provide evidence for a long-range G-type antiferromagnetic ordering commensurate with the aver¬ 
age Imma structure. Based on a general symmetry consideration, we show that the direction of 
the spins is controlled by the antisymmetric exchange imposed by the two primary structural dis¬ 
tortions, namely oxygen octahedral tilting and incommensurate atomic displacements. The tilting 
is responsible for the onset of a weak ferromagnetism, observed in magnetization measurements, 
whereas the incommensurate displacive mode is dictated by the symmetry to couple a spin-density 
wave. The obtained results demonstrate that antisymmetric exchange is the dominant anisotropic 
interaction in Fe®''"-based distorted perovskites with a nearly quenched orbital degree of freedom. 

PACS numbers: 75.25.-j 


I. INTRODUCTION 

Perovskite materials derived from the well-known 
multiferroic BiFeOa by various substitutions exhibit a 
variety of structural phases with interesting properties 
and improved functionalityii^— The electronic degree 
of freedom related to the lone pair nature of Bi^+ 
usually results in polar/antipolar atomic displacements 
in these materials. The displacements are often coupled 
to oxygen octahedral tilting. Both types of distortions 
define the symmetry of the perovskite lattice and 
control electric and magnetic properties as well as a 
cross-coupling between them. The energy landscape 
of some Bi-containing compositions consists of several 
almost degenerate phase states that can be switched 
by relatively small perturbationsJ^^— This offers a 
unique opportunity to study the structure-properties 
relationship using distinct structural modifications of 
the same material. It has been recently shown that the 
metastable perovskite BiFeo. 5 Sco. 5 O 3 can be stabilized 
in two different polymorphs via an irreversible behavior 
under heating/cooling thermal cyclingJ^ As-prepared 
BiFeo. 5 Sco. 5 O 3 ceramics obtained by quenching under 
high pressure were characterized by a complex an- 
tipolar structure with the Pnma symmetry and the 
'/2ap X dop X 2'/2ap type superstructure (op is the pseu- 
docubic unit cell). Hereafter we use p-subscript to denote 
the pseudocubic setting. The ground state atomic config¬ 
uration, however, could not be deduced unambiguously 


based on the experimental data due to existence of two 
non-equivalent Pnma isotropy subgroups different by 
the origin choice and indistinguishable in the refinement 
procedure. This configuration was determined theoreti¬ 
cally using state of the art density functional algorithms 
for structure relaxation»i^ On heating, the antipolar 
Pnma phase of BiFeo. 5 Sco. 5 O 3 was found to transform 
into the polar R3c phase identical to that of BiFe 03 . 
Subsequent cooling below the transition temperature 
resulted in onset of a novel polar phase with the Ima2 
symmetry, where the ferroelectric-like displacements of 
Bi^+ cations along the [ 110 ]p pseudocubic direction are 
combined with the antiphase octahedral tilting about 
the polar axisJi It has also been shown that both the 
Pnma and the Ima2 polymorphs of BiFeo. 5 Sco. 5 O 3 
exhibit a long-range antiferromagnetic ordering with a 
weak ferromagnetic component below about 220 K. The 
antiferromagnetic configuration was found to be of a 
G-type (where the Fe/Sc nearest neighbors in all three 
directions have antiparallel spins). 

The antipolar Pnma modification of BiFeo. 5 Sco. 5 O 3 is 
isostructural to one of the phases of La-doped BiFeOs, 
as reported by Rusakov et al. [l^ In the Bii_a;La 2 :Fe 03 
system, the Pnma phase has been found to be stable 
in a narrow compositional range and an increase of the 
La content above x ~ 0.19 induced an incommensurate 
modulation. Based on these observations, one can 
expect that the antipolar modification of BiFeo. 5 Sco. 5 O 3 
can also be driven to the incommensurately modulated 
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FIG. 1. (Color online) Electron diffraction pattern of the in- 
commensurately modulated phase of Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 . 
Indexation of the fundamental spots is done in both pseudocu- 
bic (white) and the average orthorhombic (red) cells. The 
enlarged part of the diffraction pattern (on the top) demon¬ 
strates non-overlapping second order satellites. 

structure by a partial substitution of Bi with La. 
Indeed, our present study has revealed signs of an 
incommensurate modulation in the perovskite phase of 
the Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 composition. In this work, 
comprehensive structural and magnetic studies of the 
incommensurate phase have been performed. The main 
goal was to explore how the structural modulation 
affects the spin ordering in the system. Data of the 
magnetization measurements and the neutron diffraction 
experiments were analyzed based on the symmetry 
considerations. We conclude that the magnetic structure 
of such a phase is a unique example, where a dominant 
commensurate antiferromagnetic component coexists 
with a macroscopic ferromagnetism and an incommen¬ 
surate spin density wave that has a propagation vector 
related to the structural modulation. The coupling 
mechanism has been suggested to be the antisymmetric 
Dzyaloshinskii-Moriya exchange which also fully defines 
the spin directions in the structure. 


II. EXPERIMENTAL SECTION 

High-purity oxides Bi 203 , La 203 , Fe 203 , and SC 2 O 3 
were used as starting reagents to prepare the compo¬ 
sitions of the Bii_a;Laa;Feo. 5 Sco .503 series. Previously 
calcined oxides were mixed in the stoichiometric ratio, 
ball-milled in acetone, dried, and pressed into pellets. 
The pellets were heated in a closed alumina crucible 


at 1140 K for 10 min and then quenched down to 
room temperature. The obtained material served as a 
precursor for the high pressure synthesis. The pressure 
was generated using an anvil press DO-138A with a press 
capacity up to 6300 kN. In order to avoid penetration 
of graphite from the tubular heater to the sample a 
protective screen of molybdenum foil was used. The 
samples were synthesized at 6 GPa and 1500-1600 K. 
The high-pressure treatment time did not exceed 5 min. 

An X-ray diffraction study of the powdered samples 
was performed using a PANalytical XPert MPD PRO 
diffractometer (Ni-filtered Cu Kc, radiation, tube power 
45 kV, 40 mA; PIXEL detector, and the exposition 
corresponded to about 2 s per step of 0 . 02 ° over the 
angular range 15° — 100 °) at room temperature. 

Electron diffraction patterns of the samples were 
recorded using a 200 kV JEOL 2200FS transmission 
electron microscope (TEM). The samples were crushed 
and milled with mortar and pestle. The obtained fine 
powder was dispersed in ethanol and deposited on a 
TEM grid. 

Neutron powder diffraction data for the 
Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 composition were collected 
at the ISIS pulsed neutron and muon facility of the 
Rutherford Appleton Laboratory (UK), on the WISH 
diffractometer located at the second target station)^ 
The sample (^25 mg) was loaded into a cylindrical 
3 mm diameter vanadium can and measured in the 
temperature range of 1.5 - 300 K (step 30 K, exposition 
time 2h) using an Oxford Instrument Cryostat. Rietveld 
refinements of the crystal and magnetic structures were 
performed using the JANA 2006 program^! against 
the data measured in detector banks at average 26* 
values of 58°, 90°, 122°, and 154°, each covering 32° 
of the scattering plane. Group-theoretical calculations 
were done using ISOTROPY^, ISODISTORT^ and 
Bilbao Crystallographic Server software (REPRESS 
and Magnetic Symmetry and Applications^^) . 

Magnetization data were measured using a super¬ 
conducting quantum interference device (SQUID) 
magnetometer (Quantum Design MPMS). 


III. RESULT AND DISCUSSION 
A. Crystal structure 

It was found from the obtained room temperature 
X-ray and neutron diffraction data that the crys¬ 
tal structure of Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 is different 
from the orthorhombic Pnma structure of undoped 
BiFeo. 5 Sco. 5 O 3 , indicating a compositionally-driven 
phase transition. The indexation procedure of the 
powder diffraction patterns appeared to be difficult 
using a reasonable size superstructure. Taking into 
account that about the same concentration of La has 
been reported by Rusakov et al. 0 to induce an 
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FIG. 2. (Color online) Crystal (a) and magnetic (b) structures of Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 . The crystal structure involves 
incommensurate atomic displacements (mainly Bi and oxygen as follows from the quantitative structure refinement) along the 
Uo-axis, combined with anti-phase octahedral tilting about this direction. Arrows indicate the incommensurate displacements 
of Bi (shown only for one Bi chain for clarity). The magnetic structure combines the primary antiferromagnetic component Gz 
commensurate with the average structure and two secondary ones, namely the commensurate weak ferromagnetic Fy and the 
incommensurate spin-density wave. These spin components are schematically shown as arrows along the Co-, bo- and Uo-axes, 
respectively. 


incommensurate phase in the Bii_a;Laa;Fe 03 system, an 
electron diffraction on Bio, 8 Lao. 2 Feo. 5 Sco .503 has been 
performed. The measurements confirmed the presence 


TABLE 1. Structural parameters of Bio. 8 Lao. 2 Feo. 5 Sco. 5 O 3 
obtained from the joint refinement of room-temperature 
X-ray and neutron diffraction data using the 
/mma(0,0,7)s00 superspace group with the basis vec¬ 
tors related to the parent cubic Pm3m structure as 
(1, —1, 0, 0), (0,0, —2, 0), (1,1,0, 0), (0, 0, 0,1) and origin 
at (0,1,1,0). Unit cell parameters Uo = 5.6879(1)A, 
ho = 7.9593(1)A, Co = 5.7166(1)A, = 7 c* = 0.534(3)c^ 
Reliability factors Rp = 2.79% and R^p = 3.25%. 


Atom 

Bi: 

X 

Al 

Bl 

y 

Al 

Bl 

Z 

Al 

Bl 

Uiso 

Bi/La 

0 

0.25 

0.5093(4) 

0.024(1) 


0.0475(7) 

0 

0 



-0.033(1) 

0 

0 


Fe 

0 

0 

0 

0.010(3) 


0 

0 

0 



0 

0 

0 


01 

0.25 

-0.0444(1) 

0.25 

0.037(1) 


0.023(1) 

0 

0 .012(2) 



-0.0043(7) 

0 

0 


02 

0 

0.25 

0.0785(4) 

0 .01(1) 


-0.0929(8) 

0 

0 



-0.005(2) 

0 

0 



^ and {i = x, z) are the Fourier coefficients of the first 
harmonic (n = 1) of the displacive modulation function: 

■ f^o) = I]r=o ■ fco ])+ 

B'^^cos{27Tn[rj i ■ fc^]), where rj i indicates the position of 
the j-th atom of the average structure in the l-th unit cell. 


of an incommensurate modulation with the propagation 
vector fcp = (a, a, 0; a ^ 0.27) with respect to the cubic 
perovskite unit cell (Fig. [T]). This propagation vector 
is the E-line of symmetry, following the ISOTROPY 
notationSf 22 i 2 ^ and the associated distortion is present 
in the antiferroelectric structures of the undoped 
BiFeo. 5800,503 as well as in the closely related Pbma 
structure of PbZr 03 (E 2 -distortions)i^ In both cases, 
the parameter a takes the commensurate value of 1/4 
and therefore these antiferroelectric structures can be 
regarded as the lock-in phases. Since the value of 1/4 
is not stimulated by symmetry, it can be changed by 
either external perturbations or changes in chemical 
composition. This important symmetry aspect has 
been recently highlighted by Tagantsev et al. in the 
study of lattice dynamics of PbZr 03 J^ These authors 
concluded that the antiferroelectric state is a ’missed’ 
incommensurate phase and that the transition to this 
state is driven by softening of a single polar lattice mode. 
Due to flexoelectric coupling, the system is expected 
to be virtually unstable against the incommensurate 
modulation as was shown by Axe et ali^ However, the 
Umklapp interaction forces the system to go directly to 
the commensurate lock-in phase, leaving the incommen¬ 
surate phase as a ’missed’ opportunityi^ 

This consideration provides a way to deduce the 
appropriate symmetry of the modulated phase of 
Bio. 8 Lao. 2 Feo, 5 Sco. 503 . Our experimental diffraction 
data indicate the presence of only R- and E-type 
superstructures which results in a few (3-1-1) superspace 
isotropy subgroups to be tested in the refinement 
procedurej^i^ A combination of the single-fc^ E 2 
incommensurate modulation with anti-phase octahe¬ 
dral tilting [R'^ commensurate distortions with the 
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FIG. 3. Magnetization as a function of temperature, mea¬ 
sured for Bio. 8 Lao. 2 Feo. 6 Sco ,503 under the magnetic field of 
H=100 Oe after cooling in this field (a). Magnetization loop 
measured at 5 K after a zero-field cooling (b). 


kp = ( 5 , 5 ,^) propagation vector) results in seven 
distinct superspace subgroups 22 i 2 ^ but only three of 
them are consistent with the orthorhombic metric of the 
pseudocubic perovskite unit cell. The quantitative joint 
refinement of the X-ray and neutron diffraction patterns 
confirmed the Imma{0,0,j)s00 superspace group to be 
the adequate one to describe the modulated structure of 
Bio, 8 Lao, 2 Feo, 5 Sco, 503 . The corresponding coupled order 
parameter takes the (0, S, —6; rj, 0, 0,0, 0,0, 0,0, 0,0,0, 0) 
direction in the reducible © E 2 representation 
space. The structure involves incommensurate atomic 
displacements (predominantly Bi and oxygen) along the 
[ 110 ]p pseudocubic direction combined with anti-phase 
octahedral tilting about this axis (Fig. HJa)). This 
structural model is identical to the ’missed’ incom¬ 
mensurate phase discussed by Tagantsev et al. [^| 


for the antiferroelectric PbZrOa. Earlier, Rusakov et 
al. [ 1 ^ used the /mma( 0 , 0 , 7 )s 00 superspace group 
to propose a structural model for the incommensurate 
phase of Bio. 75 Lao, 25 Fe 03 . The authors split the Bi 
position and applied a step-like occupational modulation 
to model the constant and correlated shifts of Bi and 
oxygen atoms along the Oo-axis of the average Imma 
structure. Hereafter we use o-subscript to denote the 
orthorombic setting (see Table U for the relation between 
the orthorhombic and cubic settings). In our refinement, 


FIG. 4. (Golor online) (a) Neutron diffraction patterns of 
Bio. 8 Lao, 2 Feo.sSco .503 at the vicinity of the strongest mag¬ 
netic peaks collected above and below Tjv. (b) Integrated 
intensity of the magnetic peaks as a function of temperature 
(error bars are smaller than the size of the symbol), (c) Ri- 
etveld refinement of the neutron diffraction data collected at 
1.5 K. The cross symbols and solid line (red) represent the 
experimental and calculated in the Im'ma'{0,0,'y)s00 mag¬ 
netic superspace group intensities, respectively, and the line 
below (blue) is the difference between them. Tick marks in¬ 
dicate the positions of Bragg peaks (green for satellites and 
black for fundamental). Inset shows an enlarged part of the 
diffraction pattern where the strongest structural satellites 
are observed. The solid red and blue lines represent inten¬ 
sities calculated for the incommensurate (7 = 0.534(3)) and 
commensurate (7 = 0.5) values of the modulation vector k^, 
respectively. 


the modulation was applied to the atomic displacements. 
In such an approach, the modulation is a characteristic 
of the displacive correlation function averaged over 
the sample volume (Table IJ). This model is hardly 
appropriate to analyze the local bond distances and 
angles but sufficient to consider the symmetry-controlled 
physical properties; in particular, a coupling between the 
orthogonal magnetic modes. It should be pointed out 
that although the model proposed by Rusakov et al. [l^ 
is more adequate to discuss the local crystal chemistry, 
it is another limited case which ignores the fact that the 
amplitude of the atomic displacements depends on the 
local fluctuations of La. 




























5 


B. Magnetic structure 

Magnetization measurements of Bio.8Lao.2Feo.5Sco.5O3 
revealed a weak ferromagnetic behavior below ^ 220 
K (Fig. [SKa)) with the value of the spontaneous mo¬ 
ment of ^ 0 . 022 ^B/Fe (Fig. [ 3 Kb)). This value as well 
as the critical temperature of the magnetic ordering are 
close to those in the polar and antipolar polymorphs 
of BiFeo. 5800.5031^ In agreement with the magnetiza¬ 
tion data, the neutron diffraction measurements indicate 
a long-range antiferromagnetic ordering below T-^ (Fig. 
| 4 Ka,b)). Note, a precise determination of the critical tem¬ 
perature from the neutron diffraction data is complicated 
by a superposition of the magnetic Bragg intensity with 
a diffuse component developing above Tjv. The diffuse 
component is probably related to a magnetic inhomo¬ 
geneity caused by local fluctuations in the Fe/Sc ratio. 
The magnetic Bragg reflections, observed below Tjv, are 
resolution limited and can be indexed using the aver¬ 
age orthorhombic Imma structure assuming = 0 
propagation vector. To refine quantitatively the mag¬ 
netic structure, we classified the magnetic modes accord¬ 
ing to the time-odd irreducible representations of the 
Imma space group. The crystal structure based on which 
the magnetic order emerges is however incommensurate. 
Therefore, one needs to decompose the primary struc¬ 
tural modulation in respect of the time-even irreducible 
representations of Imma. The relevant analysis revealed 
that the incommensurate structural modulation has the 
symmetry of the A4 representation associated with the 
= ( 0 , 0 ,7) line of symmetryi^ Then, by combining the 
time-odd = 0 representations with the A4 incom¬ 
mensurate order parameter, we deduced the magnetic 
superspace groups^^ and tested them in the refinement 
versus the neutron diffraction data. The toFJ © A4 com¬ 
bination resulting in the lm'ma'{ 0 , 0,7)s00 magnetic su¬ 
perspace group was found to provide a uniquely good re¬ 
finement quality (Fig. SKc)). Remarkably, this magnetic 
superspace group implies a coexistence of both commen¬ 
surate and incommensurate components of the magnetic 
order parameter (Fig. HKb)). The former is represented 
by ferromagnetic {F) and G-type antiferromagnetic com¬ 
ponents along the bo- and Co-axis of the average Imma 
structure, respectively. The latter is along the Oo-axis, 
with the modulation related to the structural one by the 
propagation vector conservation law as will be discussed 
below. The refinement procedure yielded the statistically 
significant value for the commensurate antiferromagnetic 
component only: 3 . 90 ( 4 )p,B per Fe. This value is some¬ 
what smaller than the expected one, 5 for the 
electronic configuration of Fe^’*' but practically identical 
to the ordered moment in Bio.gLao.iFeOsi^ The ferro¬ 
magnetic component in Bio.8Lao.2Feo.5Sco.5O3 is clearly 
evidenced by the magnetization data (Fig. | 3 Kb)) with 
the value well beyond the capabilities of the unpolarised 
neutron powder diffraction experiment. Although the 
presence of the modulated component does not follow 
directly from the neutron and magnetization data, it can 


be shown that the system couples it to gain full advan¬ 
tage of the antisymmetric exchange interactions. These 
interactions are anisotropic and force the direction of the 
primary G-type antiferromagnetic component to be along 
the Co-axis. 

The Heisenberg symmetric exchange interactions are 
degenerate in respect of the global spin rotations and 
therefore the corresponding part of the exchange en¬ 
ergy does not depend on the crystallographic direc¬ 
tion of the interacting spins. These interactions usu¬ 
ally dominate and, in a first approximation, define the 
relative orientations of spins in a magnetic structure. 
Then, the higher order anisotropic terms in the magnetic 
Hamiltonian such as antisymmetric exchange, single ion 
anisotropy, magnetoelastic coupling and dipole-dipole in¬ 
teractions remove the global spin rotation degeneracy, 
making some crystallographic directions preferable. In 
the Bio.8Lao.2Feo.5Sco.5O3 perovskite, the octahedrally 
coordinated Fe^’*' ions with a nearly quenched orbital de¬ 
gree of freedom [L — 0,3 = 5 / 2 ) interact antiferromag- 
netically via the strong superexchange, which results in 
the G-type magnetic structure as found experimentally 
and as observed in many other Fe^“''-based perovskitesi^ 
Structural distortions activate the anisotropic terms; in 
particular, antisymmetric exchange. 

As follows from the structure refinement, described in 
the previous section, there are two primary structural 
distortions in Bio.8Lao.2Feo.5Sco.5O3, namely anti-phase 
octahedral tilting and incommensurately modulated Bi 
and oxygen displacements. The part of the antisymmet¬ 
ric exchange related to octahedral tilting was considered 
in details in Ref. [i3- It has been shown that the ax¬ 
ial distortions associated with the octahedral rotations 
are responsible for the weak ferromagnetic properties of 
antiferromagnetically ordered perovskites with a G-type 
spin configuration. The relevant part of the Dzyaloshin- 
skii vector, is expressed by the antiferroaxial vec¬ 

tor which is a characteristic of the tilting pattern. In 
other words, the component of the Dzyaloshinskii vector 
which induces the weak ferromagnetism coincides with 
the tilting axis of octahedra. One can therefore expect 
that the spin components of the primary G-type antifer¬ 
romagnetic mode are confined within the (boCo) plane, 
since the octahedra are tilted about the Oo-axis in the 
average Imma structure. This part of the antisymmetric 
exchange energy, D°j ■ [Si x Sj], is degenerate in respect 
of the moment direction in the {bgCo) plane. This degen¬ 
eracy is however removed, when we take into account the 
antisymmetric exchange imposed by the incommensurate 
structural distortion. 

To demonstrate that, we need to work out the appro¬ 
priate free-energy coupling terms in respect of the parent 
Pm 3 m symmetry (see Appendix for details). The incom¬ 
mensurate structural distortion transforms as the twelve¬ 
dimensional Yi2{r]j=i-i2) irreducible representation asso¬ 
ciated with the fcp = (a, a, 0) propagation vector of the 
cubic space group (a = ^). The direction of the or¬ 
der parameter in the Eg representation space is specified 
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by the single non-zero component rji. The symmetry 
properties of the G-type antiferromagnetic mode are de¬ 
fined by the = (^, 5 , 5 ) propagation vector and the 
three-dimensional representation A com¬ 

bination of these order parameters to a free-energy invari¬ 
ant requires a coupling to a time-odd (magnetic) physical 
quantity, with the modulation related to the structural 
one as fc™"® = (^ — a, ^ — a, i) to maintain the transla¬ 
tional invariance. As detailed in the Appendix, a third 
power invariant with a magnetic dipole order parameter 
can be formed only when the spins in the primary G-type 
antiferromagnetic mode are along the orthorhombic Cq- 
axis (/ii = 0 ,/i 2 = /Ta 0). The relevant energy term 
is: 


M2??iCio + = 2/r77^ (1) 


which describes a coupling of the spin density wave with 
the spin components being along the Oo-axis of the av¬ 
erage Imma structure (Fig. HKb)). The presence of this 
incommensurately modulated spin component is in full 
agreement with the magnetic Im'ma' super¬ 
space symmetry derived above. In the framework of the 
representation theory, the symmetry properties of the 
spin density wave are specified by the twelve-dimensional 
time-odd irreducible representation with 

the single non-zero component of the order parameter 
^ 10 . The crucial point is that there is no such coupling, 
if the spins in the antiferromagnetic configuration are 
along the 60 -axis (/j,i ^ 0 ,/i2 = Ms = 0). Thus, the 
system chooses the Cg-axis for the spins direction in 
the primary antiferromagnetic G-mode to activate the 
energy term specified by expression o, which breaks 
the degeneracy between the ho- and Co-axes. 

The free-energy term specified by expression m is not 
invariant under a global spin rotation (it vanishes when 
the spins are along the orthorhombic 60 -axis) revealing 
its relativistic nature^i This invariant implies that 
the incommensurate atomic displacements modulate 
the relevant component of the Dzyaloshinskii vector 
which in turn induces the spin density wave through the 
relativistic antisymmetric exchange. 

Thus, the magnetic structure of Bio. 8 Lao. 2 Feo. 5 Sco ,503 
can be fully understood by taking into account only 
the isotropic symmetric and anisotropic antisymmetric 
exchange interactions. The first type of interactions 
defines the primary G-type antiferromagnetic config¬ 
uration through the strong 180-degree superexchange 
expected for the half-occupied eg orbitals of Fe^’*'. The 
second (anisotropic) type of interactions chooses the 
Co-axis for the spins direction in the primary mode to 
fully exploit the dominant structural distortions and 
couple the secondary orthogonal ferromagnetic and 
incommensurate spin components. 


IV. CONCLUSIONS 

A 20% substitution of Bi with La in the 
Bii_ 3 ;Laa;Feo. 5800,503 system synthetized under 
high-pressure and high-temperature conditions in¬ 
duces incommensurately modulated structural phase. 
The symmetry of this phase is described by the 
/mma(0,0, 7 )s 00 superspace group (7 = 0.534(3)) with 
modulated displacements of Bi/La and oxygen ions. The 
structure combines the same type of primary distortions, 
as the prototype antiferroelectric structure of PbZrOs. 
The difference between the two structures is in the 
value of the propagation vector for the antipolar dis¬ 
placements (E 2 distortive mode) which is commensurate 
(7 = 0.5) in the case of PbZrOs. Both propagation 
vectors (commensurate and incommensurate), however, 
belong to the same line of symmetry and therefore the 
commensurate value (lock-in phase) is not stimulated 
by the symmetry and can be tuned by composition 
or external perturbations such as pressure, strain (for 
thin-film forms) and electric field. 

Below ~ 220 K, a long range antiferromag¬ 
netic ordering commensurate with the average Imma 
structure takes place. The spins are aligned along the 
Co-axis which allows the system to gain energy from 
the antisymmetric exchange activated by the two pri¬ 
mary structural distortions: namely, octahedral tilting 
and the incommensurate atomic displacements. The 
antisymmetric exchange imposed by the tilting induces 
a weak ferromagnetic component along the 60 -axis. 
The part of the antisymmetric exchange related to the 
modulated atomic displacements with = {a, a, 0 ) 
couples a spin-density wave with the propagation vector 
= (| — a, ^ — a, 5 ) and the spin components 
being along the Uo-axis. These results demonstrate 
the crucial role of the antisymmetric exchange in mag¬ 
netic properties of Fe^“'"-containing distorted perovskites. 
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APPENDIX 

To explore the Landau free-energy terms activated by 
the incommensurate structural modulation, let us de¬ 
fine the transformational properties of the distortions in¬ 
volved in respect of the parent Pmim symmetry. In this 
cubic structure, Bi/La, Fe/Sc and oxygen atoms occupy 
16(i, i, i), la(0, 0, 0) and 3d(i,0,0) Wyckoff positions, 
respectively. The structural modulation and the G-type 
spin configuration transform as the twelve-dimensional 
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time-even E 2 (r?i, r? 2 , %, ??4, ??5, ?76, ?77, ??8, %, ?yio, ?yi i, ?? 12 ), 
{fcp = (a, a, 0)} and the three-dimensional time-odd 

/i 2 , Ms), = (^> 5 , 5 )} irreducible repre¬ 

sentations of the cubic space group ({} indicates a wave 
vector star with a representative arm enclosed)The 
lowest degree coupling term maintaining the time re¬ 
versal symmetry is a cubic trilinear invariant of the 
j I ~ type- Here, ^(-denotes the coupled time- 

odd order parameter whose symmetry we need to fig¬ 
ure out for the cases when the magnetic moments are 
along the orthorhombic bo (pseudocubic [OOlJp) and Co 
(pseudocubic [110]p) axes. These spin configurations oc¬ 
cur when the mRl order parameter takes the (/ri, 0 , 0 ) 
and {0, ^ 2 , fJ '2 = Ms) directions, respectively. The 
incommensurate atomic displacements along the [ 110 ]p 
pseudocubic direction (orthorhombic Oo-axis) are de¬ 
scribed by the (ryi, 0,0, 0,0, 0,0, 0,0, 0,0, 0) order pa¬ 
rameter in the S 2 representation space. The trans¬ 
lation symmetry requires the ^-quantity to be associ¬ 
ated with the = (| — a,^ — Q!,i)} propagation 

vector star (S'-line of symmetry), where a = 7/2 is 
the wave number of the structural modulation. This 
comes directly from the trilinear form of the invariant 
which requires the product of the Fourier transforms 
g-27rz(fep -t) g27i-*(fej, -t) associated with the ijj and 

order parameters to change sign at the tpi — ( 1 , 0 , 0 ), 
tp 2 = (1,0,0) and tp^ = (1,0,0) translations. There 
are four twelve-dimensional irreducible representations, 
m5'i,=i_4(^l, ^ 2 , ^31^4, ^51^6; ^7; Cqi ^ 10 ; Cll; C 12 ): aSSO- 
ciated with the star, but only three of them ap¬ 

pear in the decomposition of the pseudovector (magnetic) 
reducible representation localized on the Fe 16 Wyckoff 
position: 


Tmag{^b) = mS2®rnS3®mS4 ( 2 ) 

Using the ISOTROPY software (irrep version 2011) , 
one can derive that the order parameter with the mSi 
symmetry is coupled to the (g) S 2 product through 

the general free-energy invariant of the form: 

AiibiCg - + 

IJ‘2r]9^5 — ^2^10^6 — M27?11^7 + ^2^12^8 + 

^J-3V5^1 - /i37?66 - f^3V7^3 + Ai3b8?4 (3) 


which is reduced down to the simple term for the 

relevant direction of the E 2 order parameter. This term 


describes the allowed coupling scheme for the case of the 
G-type spin configuration with the moments being along 
the orthorhombic 60 -axis (pseudocubic [001]p). The mSi 
irreducible representation, however, does not enter into 
the decomposing of the reducible pseudovector represen¬ 
tation on the Fe Wyckoff position (zero subduction fre¬ 
quency, see expression ©)• This means that there are no 
any dipole magnetic modes localized on the Fe site with 
the symmetry of the mSi representation and therefore 
the system cannot activate the antisymmetric exchange 
through this type of the tri-linear invariant. 

The situation is different when the coupled order pa¬ 
rameter transforms as the mS 3 irreducible representa¬ 
tion. The general coupling invariant takes the form: 


^J'lV5^2 + MWeCi - - ^J-lV8^3 + 

+ MlbloCs — Mlbll ?8 — Mlbl2^7 + 

M2bl?10 + M27?2^9 + /427?3 Ci 2 + M2b4Cll + 
f^2V5^2 + M 2 b 6 Cl + ^2777^4 + ^J-2VS^3 + 

M3?7l?10 + M37?2^9 — M3b3?12 — M3b4Cll + 

^ J - 3 V 9^6 + ^37710^5 + M377 i 1^8 + M37712C7 (4) 

with the non-vanishing terms ^ 2111^10 + MsbiCio for the 
(771,0,0,0,0,0,0,0,0,0,0,0) direction describing the in¬ 
commensurate structural modulation. These terms are 
activated when the spins of the primary G-type antifer¬ 
romagnetic configuration are along the orthorhombic Cq- 
axis (pseudocubic [110]p). The subduction frequency of 
the mSs representation is non-zero for the 16 pseudovec¬ 
tor reducible representation ex. @ and the projection op¬ 
erator yields the spin density wave localized on the Fe 
site with the spin direction along the orthorhombic Oq- 
axis (pseudocubic [110]p) (Fig. Oi). 

Finally, let us point out, for completeness, that the 
magnetic order parameter transforming as the mS 2 irre¬ 
ducible representation does not form a third power invari¬ 
ant with the mR^ 0 E 2 product. The mS '4 representation 
does form such invariant and provides coupling, when the 
spins, in the primary G-type configuration, are along the 
orthorhombic Co-axis. The latter situation is however 
unfavourable for the part of the antisymmetric exchange 
associated with the octahedral tilting as discussed in the 
main text. 
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